Slow quench dynamics of periodically driven quantum gases 
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We study the evolution of bosons in a periodically driven optical lattice during a slow change of the 
driving amplitude. Both the regime of high frequency and low frequency driving are investigated. In 
the low frequency regime, resonant absorption of energy is observed. In the high frequency regime, 
the dynamics is compared to a system with an effective Hamiltonian in which the atoms are 'dressed' 
by the driving field. This 'dressing' can dramatically change the amplitude and sign of the effective 
tunneling. A particular focus of this study is the investigation of the time-scales necessary for the 
evolving quantum state to follow almost adiabatically to the ground-state of the effective many body 
system. 

PACS numbers: 03.75.Kk, 05.70.Ln, 03.75.Lm, 67.85.Hj, 
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I. INTRODUCTION 

In experiments with cold atomic gases confined in op- 
tical lattices, many parameters can be tuned almost at 
will [lj. This makes them a very promising set-up in or- 
der to study long standing questions as the origin of high 
temperature superconductivity, or the existence of spin 
liquids in frustrated quantum magnets. 

Recently the use of high frequency periodical driving 
forces has been proposed to engineer complex effective 
many body Hamiltonians not reachable by other means. 
The effect of periodic forces to 'dress' atoms is well known 
from quantum optics where potentials for atoms can be 
designed by the application of light fields 0, Q. The 
dressing is not only efficient for non-interacting atoms. 
For weakly interacting Bose-Einstein condensates, an ef- 
fective suppression of tunneling in double well structures 
4] and also in one-dimensional optical lattice potentials 
5| as been observed. In the context of strongly inter- 
acting atom gases, interesting effects of high frequency 
driving forces are predicted as for example the change 
of sign of the tunneling amplitude in an optical lattice 
(| 0] which can induce an effective negative temperature 
8(, controlled bound-pair transport [9(, the building of 
frustrated antiferromagnetism at experimentally feasible 
temperatures [10( and the generation of effective stag- 
gered magnetic fields and hence unconventional super- 
conductivity [ll| . An important step towards the design 
of such complex quantum phases by driving potentials in 
strongly interacting systems was the recent experimental 
realization of the transition from a superfluid to a Mott- 
insulator by accelerating the lattice [l2]. This proof of 
principle experiment proposed in Ref. [6[ shows the fea- 
sibility of the generation of effective models in atomic 
gases with strong interactions despite the occurrence of 
many-body resonances. 

However, important questions to address are not only 
the engineering of complex effective Hamiltonians, but 
also how well for example the ground state of the effec- 
tive Hamiltonian could be reached by slowly switching on 
the driving potential and how much heating the driving 



potential induces. To answer these questions the almost 
perfect isolation of cold atom gases by which no thermal- 
ization with any surrounding environment takes place has 
to be taken into account. Due to this isolation the prepa- 
ration of the atomic cloud is almost entropy conserv- 
ing (and temperature is not conserved), if the changes 
of parameters are performed slowly enough. The ques- 
tion of whether it is possible to follow the ground-state 
of a strongly correlated state while slowly changing the 
parameters of the non-periodically driven system is still 
an unresolved puzzle [13l . Il4j . In the context of bosonic 
atoms confined in an optical lattice, it has been addressed 
recently both theoretically |15l — TlTl j and experimentally 
\m\m- The presence of many avoided and unavoided 
level crossings in the many-body spectrum makes an un- 
derstanding very involved. For the crossing of a phase 
transition to an ordered phase the dynamics has been 
proposed to obey a universal mechanism, the so called 
Kibble- Zurek mechanism [20l. l2lj]. This mechanism pre- 
dicts an algebraic dependence of the defect generation on 
the ramp time. However, its validity in quantum systems 
is not yet established [la . [lij . 

Even more involved is the question of adiabaticity in a 
periodically driven many body system [22;]. The descrip- 
tion of the periodic driving is often performed using the 
so-called Floquet formalism [2^, [24| . This formalism uses 
the fact that due to the periodicity in time, a Brouillon- 
zone-like structure occurs in the quasi-energy spectrum. 
Whereas in a non-driven system the energy values follow 
smoothly a change in a time- varying parameter, the Flo- 
quet quasi-energies are found to exhibit discontinuities in 
the thermodynamic limit during the change of the ampli- 
tude of the driving [22j. Therefore strictly speaking an 
adiabatic limit does not exist. However for many physical 
relevant situations of large system sizes the quasi-energies 
are expected to remain close to the values of the undriven 
system [22| . 

Previous numerical studies in many body systems have 
investigated the behavior of a periodically driven system 
for small system sizes (approx. 12 sites) Q. Different ex- 
citations depending on the driving frequency have been 
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identified due to the crossing of (avoided) level cross- 
ings in the energy spectrum. In the regime of high fre- 
quency driving, a decent following of the ground state 
value has been pointed out. An improvement of the adi- 
abatic following has been suggested by a shortening of 
the turn-on time, since in this way, at narrow crossings, 
Landau-Zener transitions could be suppressed |6j. 

It is not clear how these findings in small systems carry 
over to larger system sizes in which finite size gaps de- 
crease and dense energy bands develop. In our work we 
will investigate the response of a bosonic atom cloud of 
realistic size confined to an optical lattice potential to a 
periodic driving while slowly increasing the driving am- 
plitude to a desired value. As a driving force we consider 
the tilting of the optical lattice potential which exper- 
imentally can be generated by the use of a detuning of 
the frequencies @ or a spatially oscillating mirror reflect- 
ing the lattice beams [25| . We investigate three different 
physical situations: (i) the tuning to a vanishing tunnel- 
ing amplitude, (ii) the decrease of the tunneling ampli- 
tude to a final value without changing the sign and (iii) 
the inversion of the sign of the tunneling amplitude. This 
choice is motivated by different physical situations. Due 
to the direct proportionality between the force induced 
and the mass of the atoms, the driving could be used as 
an efficient method to manipulate the tunneling constant 
of only selected species in an optical lattice (i.e. the heav- 
ier ones). Thereby generating a vanishing effective tun- 
neling (i) could be used to create a quenched disordered 
potential for the heavier species. Situation (ii) is a more 
general change of a parameter in a system, interesting, for 
example, in the study of a superfluid to Mott-insulator 
transition. Finally situation (iii), i.e. the inversion of the 
sign of the tunneling amplitude, is a very interesting tool 
to engineer more complex quantum phases as spin liquids 
in frustrated geometries Our study of the inversion 
of sign in a simple one-component bosonic system will be 
the test-bed if this inversion is possible. In particular we 
investigate whether it is possible to produce a state which 
is close to the ground-state of the effective Hamiltonian 
without producing too many excitations when inverting 
the sign of the tunneling parameter. 



II. SETUP 

In experiments with cold atomic gases, potentials of 
different geometries can be engineered by the application 
of light or magnetic fields. Here we consider the applica- 
tion of two far-detuned counter-propagating laser beams 
to a one-dimensional tube of bosonic atoms. The atoms 
feel a potential V(x,t) of the form: 



V(x,t) = Vbsin 2 (k L (x + X(t))) . 



(1) 



Here fcz, = 2it/\l is the wave vector and Xl the wave- 
length of the laser, t denotes the time. The potential 
strength Vq can be tuned experimentally, since it is pro- 
portional to the intensity of the laser. Additionally to 



the standing wave usually created to form a steady peri- 
odic potential, k^Xl!) is a time-dependent phase which 
can be adjusted by introducing a dephasing of the two 
counter-propagating lasers, e.g. by a moving mirror. 
We consider that the phase varies periodically in time, 
i.e. X(t) = A(t)cos(uit) (where A(t) is the amplitude of 
the oscillations). 

Considering a sufficient height of the periodic poten- 
tial, many physical properties of the ultracold bosonic 
atoms can be described within a one-band Bose-Hubbard 
Hamiltonian @, [H, 13 



H — H tun + H int + (t) ^ 



a j rij 



(2) 



where H+- 



h.c. 



U 



and H int = — ^ - l ) 



The operator bj (St) annihilates (creates) a boson at 

site j. The operator hj = by>j is the number operator on 
site j and a denotes the lattice spacing. The first term 
models the kinetic energy of the atoms with the tun- 
neling amplitude J and the second term the short range 
interaction between the atoms with the onsite interaction 
strength U . The third term, which describes a force peri- 
odically varying in time with frequency uj, stems from the 
detuning of the counter-propagating laser beams. This 
force is related to the phase oscillations of the detuned 
laser beams by 



F u (t) = -rn 



d 2 X(t) 
dt 2 



(3) 
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ui 2 A(t) cos(wi) 
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where too is the mass of an atom. If the amplitude A(t) 
is varying slowly compared to uj the expression takes the 
simpler form [28[ 

F u (t) = /(i)cos(urf) with f(t) = A(t)m uj 2 . (4) 

We will investigate this model for realistic system sizes L 
around 48 to 64 sites with open boundary conditions. We 
use the time-dependent densit y m atrix renormalization 
group method (t-DMRG) [29l433j keeping up to a few 
hundred states if necessary. 

Before we turn to the description of the driven Bose- 
Hubbard model, let us summarize some of the proper- 
ties of the Bose-Hubbard model in the absence of the 
driving, i.e. F u = 0. For weak interaction strength 
(U/J < (U/J) c ) a superfluid state is present, in which 
atoms are delocalized over the system. In one-dimension, 
two main characteristics of this state are an asymptotic 
algebraic decay of the single-particle correlation func- 
tions with distance and a gapless excitation spectrum. In 
contrast at strong interactions (U/J > (U/J) c ) a Mott- 
insulating state occurs at commensurate filling. In this 
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state the atoms are localized in space, i.e. the single- 
particle correlation functions decay exponentially with 
distance and particle number fluctuations are strongly 
suppressed. The Mott-insulating state is incompressible 
and has a gap to the first excitations. For large inter- 
action strengths, the gap is of the order of the interac- 
tion amplitude U. At incommensurate filling the system 
exhibits a crossover between the superfluid and the so- 
called Tonks-Giradeau gas at large interaction strengths. 
Note that in this limit of strong but finite interaction 
strength, many low energy properties can be described 
by a gapless fermionic system. However, in dynamical 
situations higher energy bands can play an important 
role. 

The approximation of taking only a single-band Bose- 
Hubbard model ([2]) is only valid as long as excitations to 
the higher lying bands can be neglected. In particular 
an energy quantum huj of the driving field must be far- 
resonant with the excitation energies to the next bands. 
In the harmonic approximation, the band gap between 
the lowest and first excited Bloch band of the lattice is 
given by ujho ~ 2Er^JVq/ -Er. Here we will only consider 
driving frequencies uj lower than LO} lo . 



III. FLOQUET DESCRIPTION AND 
EFFECTIVE HAMILTONIAN 

Neglecting the slow increase of the driving amplitude, 
the considered system falls into the wide class of periodic 
Hamiltonians H(t) — H(t + T), where T is the period 
and u) = 2w/T the frequency. This periodicity leads to 
specific properties which can be discussed within the so- 
called Floquet formalism [23|, [24[ . In the case of a fast 
driving (we will specify later on what we mean by fast) an 
effective time-independent Hamiltonian can give a very 
good description of many physical properties. In the con- 
sidered case at a fixed time t this effective Hamiltonian 
Q is given by 

Heft = -JJ (f(t)a/hw) ^(StS i+1 + h.c.) 

3 

3 

where J n denotes the Bessel function of n-th order. Com- 
paring to the Hamiltonian @ without driving an effective 
tunneling amplitude has been generated by the driving 
which we will denote by J c s(t) = J Jo{f{t)a/hoj). 

We recapitulate in the following the derivation of this 
effective Hamiltonian, since it gives many useful insights 
into the behavior of the system also outside the range of 
validity of the effective model. 

For a Hamiltonian periodic in time, a complete set of 
solutions to the Schrodinger equation (h - ihd t ) W) = o 
is given by the states \ip) = \u n {t))e~ lEnt / h where \u n (t)) 



are T— periodic solutions of the eigenvalue problem 

fl(t)Mt)) = E n \u n {t)). (6) 

Here $j(t) is the Floquet Hamiltonian which is given by 
Sj(t) = H — ihdt, E n are different quasi-energies and 
\u n (t)) are the so-called Floquet states. If \u n (t)) is a 
Floquet state, the state \u n (t),m) — \u n (t))e muit (with 
m integer and u) = 2tt/T) is also a T— periodic solution 
of equation ((6]) with quasi-energy E n + mhw. Hence, in 
analogy to spatially periodic systems, the quasi-energy 
spectrum has a repeating Brillouin-like structure with 
the width of the Brillouin zone being tko. 

The aim is to find a time-independent effective Hamil- 
tonian, H c ff, with an energy spectrum approximating 
the quasi-energy band E n G [— hoj/2, hcj/2] of the Flo- 
quet Hamiltonian Sy (t) averaged over one period of time. 
In order to derive this effective Hamiltonian, one per- 
forms a transformation of the basis states \u n (t)) to a 
new opportunely chosen basis. In general, the new ba- 
sis can be built using a T— periodic Hermitian operator 
F(t) = Ft(t) = F(t + T) by \n(t),m) = U F , n {t)\u n {t)) 

with U F , m {t) = e - lp ^ +lmut and m G Z. The matrix 
elements of the Floquet Hamiltonian time-averaged over 
one period become |34| : 

((n(t),m\Sj(t)\n'(t),m')) T = 

6 mim ,(n\(mhuj + (S(*) (0,0) )tK} (7) 

+ (1 - ^, m ')(«l(e l(m '" mM i3W^ 0) > T |n') (8) 

Here we used the notation f)(i)( m,m ) = 

U F m (t)F)(t)UF,m>(t) and (...}t denotes an average 
over one period T. The energy spectrum of the 
Hamiltonian H c ff = (Jo(t) < - 0,0 - l }T well approximates the 
quasi-energy spectrum of the time-averaged Floquet 
Hamiltonian f) (t) if the coupling between different bands 
can be neglected. 

In the considered situation Eq. ©, a good choice 
of the operator is F(t) = (fa/hu) sin(wt) J^ - jrij. 
Starting from the basis formed by the Fock 
states the new basis becomes |{n;},m) = 

|{ ni }) e -«(/a/^)sm( W t)E J in i +im W t_ The matl]x elements 

of the time-averaged Floquet Hamiltonian become 

(({n' l },m'm)\{ni},m)) T = 

S m ' ,m{{n'i\\mfru) + H- mt + joH tun \{ni}) (9) 

+ (1 ~ 5 m \ m )3s{m-m'){{n'i}\H tun \{ni}) (10) 

where j a = J a (fa/hoj) and s = £\ j(rij -n'j). 

We see that if the non-diagonal coupling vanishes, 
Hu> 3> U or J, and only one zone E n E [—hu)/2,hoj/2] 
is taken into account, we recover the effective Hamilto- 
nian ([5]). This means that this effective Hamiltonian de- 
scribes well the properties of the system if no absorption 
event described by the off-diagonal part takes place. It 
is indeed important to understand, in experimentally in- 
teresting problems for one-dimensional gases, when this 
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approximation is valid. We will discuss the non-diagonal 
term more in detail in the next section, since beside its 
importance for the validity of the effective Hamiltonian 
it also reveils a lot of information on the resonant absorp- 
tion at low frequencies. 



IV. RESONANT VERSUS EFFECTIVE 
BEHAVIOR 

In this section we discuss the behavior of the system 
subject to a driving with a wide range of frequencies. Mo- 
tivated by recent experimental work 0, [l2| , in which the 
driving amplitude was increased linearly with time, we 
will -if not stated otherwise- focus on a linear increase 
of the amplitude of the phase detuning f(t) = f t t. There- 
fore we do not deal with an exact time-periodic Hamil- 
tonian anymore. However, we found that if the change 
of the amplitude is very slow compared to the driving 
frequencies, the properties of periodic systems are still a 
reasonable approximation. 

A typical response of the system subject to the driv- 
ing, is summarized in Fig. [TJ Note, that these are results 
of the exact time-evolution using the t-DMRG method. 
Fig. [T] shows how the compressibility Kj = (n?) — (fij) 2 
at j = L/2 evolves in time when the driving amplitude 
increases. In order to have a lighter notation, we will use 
k = kl/2- The final driving amplitude is chosen, such 
that the final value of the effective tunneling vanishes 
(case (i)). At the final time of the evolution the com- 
pressibility of the ground-state of the effective Hamilto- 
nian (from now on we will refer to the ground-state of 
the effective Hamiltonian at the final driving strength 
as the target-state) is zero. Let us also note that the 
compressibility at uniform filling is related to the inter- 
action energy of the system by Ej = U/2 J2j K j + const. 
Therefore the absorption of interaction energy is directly 
reflected in the compressibility. 

For low driving frequencies a resonant absorption of 
energy, i.e. a strong increase of k, in certain frequency 
bands occurs (cf. Fig. [2] hut/ J = 6,14). In contrast for 
large frequencies a slow decrease of n is visible. This slow 
decrease appears to be almost independent of the driving 
frequency (cf. Fig. |2J huo/J = 18 and hw/J = 30). This 
is the regime in which the behavior of the system can be 
described by the effective model (JS|) taking at each time 
t the effective tunneling amplitude J c s(t). The effective 
tunneling amplitude slowly decreases with time causing 
a decrease in the compressibility. This effective model 
compares reasonably well to the exact numerical calcu- 
lations at high frequency (cf. Fig. Let us note that 
for some frequencies the time-evolution shows a mixed 
behavior between the two discussed cases. A crossover in 
time might occur between the decrease of the compress- 
ibility for small time, i.e. the diagonal term dominates, 
and the absorption of energy for larger time, i.e. the non- 
diagonal term dominates (cf. Fig. [21 haj/J = 12). 

In the following we will discuss the occurrence of the 
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FIG. 1. (color online) Time-evolution of the compressibility k 
for different driving frequencies ui. The value of the strength 
of the driving force f(t) is linearly increased in time, such 
that at time tf = 12h/J, ^(fttfa/hui) = 0. The 'undressed' 
interaction strength is U/J = 4 and system length L — 48 
with N = 48 atoms. 
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FIG. 2. (color online) Time-evolution of the compressibil- 
ity k for chosen driving frequencies w (ftw = 6J, stars; 
liuj — 12J, squares; fiw — 14J, triangles; liui — 18J, circles; 
liw = 30, diamonds) and without driving following the de- 
crease of the effective tunneling amplitude (continuous line). 
The inset shows the evolution of the effective tunneling am- 
plitude. Same parameters as in Fig. [1] 



resonances. Similar resonances have already been studied 
in different contexts [35T - [37l | . In the considered model 
in the strong coupling limit resonant behavior has been 
investigated by linear response theory [38j and in the 
Floquet formalism [?J investigating the off-diagonal term 
(Eq. (fTU)Q in perturbation theory. 

A strong resonance occurs at U « fko corresponding 
in the large U limit to a single particle-hole excitation. 
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Weaker resonances describing multi-photon processes to 
the same single particle-hole excitation are expected at 
mtvjj « U (m is an integer). The resonant behavior 
around U ~ huj is clearly seen in Fig. [TJ However, due 
to the chosen intermediate value of U j J = A > (U/J) c 
this resonance is very broad [39| and no higher photon 
resonance can be resolved. The coupling to the naively 
expected two particle-hole excitations at mfku 2U van- 
ishes up to second order in J/U for systems at filling one 
0, HHdl]. In contrast, excitations at 3U are expected. 
In second order in J/U, only resonances at mftw ~ ?>U , 
where m is an even integer, have non- vanishing matrix 
elements 0. In Fig. [T] a resonance at 2Kuj ~ 3U seems 
to be present leading to a substructure of the broad res- 
onant region. Additionally to the predicted resonances, 
in Fig. [TJ a weaker resonance is seen which shifts from 
huj ~ AU towards huj w 3U. In the strong coupling limit 
this resonance corresponds to a process in which two par- 
ticles from a singly occupied site hop to the same site of 
filling one between them (thus forming a triply occupied 
site). A shift of this resonance with the parameter U/J 
has already been found in linear response calculations in 
the context of the modulation spectroscopy [37[ ■ Hence 
the resonance is related to the spectral properties of the 
undriven system. We further find that a weak branch sep- 
arates from this weaker resonance and joins the frequency 
band around hui = 2U. We believe that this branch is 
due to subtle effects of the dressing of the original energy 
states by the strong driving. Even though resonances at 
higher value of the frequency Huj exist if higher orders in 
perturbation theory are considered, they seem to be very 
much suppressed for hui > bU at the shown parameters. 
In this regime the off-diagonal coupling terms (|10j) can be 
neglected and the effective description approximates the 
behavior of the system. The range of frequencies above 
which the effective regime is valid depends very much 
on the system parameters. For example at unit filling 
and low interaction strength U/J « 1 — 2 the frequency 
range above hui = 12 J is well described whereas at larger 
interaction strength, Huj > 5U is a relatively safe approx- 
imation. 

The behavior at resonance changes for different values 
of the filling [36| . We found additional resonances and 
shifts in their position and weights at larger fillings. For 
example at filling n — 1.5 and n = 2, a strong resonance 
occurs in the compressibility at Huj ss 5U. In the strong 
coupling limit this resonance can be associated with the 
tunneling of two atoms in a single occupied site to a dou- 
bly occupied site located between them. This forms a 
site with four atoms with its two neighboring sites non- 
occupied and has an energy difference to the ground-state 
of about 5U. This process is highly suppressed for unit 
filling but becomes important for larger fillings since the 
probability to have a doubly occupied site strongly in- 
creases. In the experiments it will be important to avoid 
these additional resonances either by detuning the driv- 
ing frequency or by assuring that the density in the trap- 
ping potential does not exceed a certain filling. 



Let us note that the resonant energy absorption de- 
stroys very quickly the coherence between the atoms 
which can be seen by a rapid decrease of the single- 
particle correlations (not shown). 



V. SLOW QUENCH DYNAMICS IN THE HIGH 
FREQUENCY REGIME 

We now focus on the regime of high frequency driving. 
We aim to quantify the quality of the preparation of a 
target-state. The situation we consider is a slow linear 
increase of the driving amplitude f(t) — ftt over a finite 
ramp time tf. This preparation by almost adiabatic in- 
crease is intriguing since in infinite systems the Floquet 
states and their quasi- energies have been shown to be 
discontinuous functions of the driving amplitude almost 
everywhere [22l |. 

Previous numerical studies focused on small system 
sizes of at most 12 sites [1, 0j- In these small systems, 
finite size gaps occur and could lead to a very different 
behavior of the dynamics. We study the dependence of 
the proximity of the final state to the target-state on 
the ramp time tf in large, and experimentally relevant, 
systems. To distinguish the effect of the driving force 
from the adiabatic change in the effective tunneling am- 
plitude, we compare the driven situation to the response 
of the system to a slow change of the tunneling ampli- 
tude without the presence of the driving (at each time 
the tunneling amplitude is chosen to correspond to the 
effective tunneling amplitude which would have been in- 
duced by the driving). The later situation will be called 
'undriven' in the following. Three physically different 
situations are distinguished (i) the preparation of a state 
with the vanishing effective tunneling (J e f{(tf) = 0), (ii) 
a state corrsponding to a reduced but finite effective 
tunneling (J e g(tf) — J/2), and (iii) a state in which 
the sign of the tunneling amplitude has been reversed 

(J eS (tf) « -0.4J). 



A. Slow decrease of the effective hopping when 
just in the Mott region 

The initial state we consider here is the ground state 
at the interaction strength U/J = 4 at the filling of one 
particle per site. This state lies just above the phase tran- 
sition in the Mott-insulating phase (the transition, for a 
homogeneous system at unit filling, is at U/J w 3.38 
[4l[ |42j ) . In Fig. [3] the final value of the compressibility 
k versus the final time of the effective quench is shown 
for the situations (i)-(iii). The behavior of different driv- 
ing frequencies uj is compared to the expected behavior 
from the undriven case. In an experiment the reduction 
of the tunneling amplitude can be directly performed by 
the increase of the optical lattice strength. However, the 
change of sign cannot be induced by the simple modula- 
tion of the lattice height. 
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FIG. 3. (color online) Compressibility k versus the final 
time of the quench tf for different driving amplitudes ft 
such that (a) J e s(tf) = , (b) J c g(tf) — 0.5J, and (c) 
J e s(tf) = — 0.4J. Different curves are hw/J — 18 (dash- 
dotted line) and hw/J = 25 (continuous line), undriven 
(dashed line). The symbols denote actual numerical data. 
The horizontal lines represent the values of the compressibil- 
ity of the target-states and arrows indicate the value of the 
compressibility in the initial state. Parameters are U/J = 4, 
N = L = 48. 



In the undriven system, for short ramp times tf the 
value of the compressibility shows a rapid decrease from 
the value of the initial state with the ramp time [43| . 
For long ramp times tf the value of the compressibility 
reached in the undriven system is in all situations close to 
that of the target-state. In the driven systems the com- 
pressibility shows roughly a similar decrease. However, 
different features have to be commented. For case (i) and 
(ii) this decrease of the compressibility agrees well with 
the undriven system for small ramp times tf, but devi- 
ates considerably for larger ramp times tf. For a very 
large frequency huj — 25 J an almost monotonic decrease 
is found for large ramp times, whereas fku — 18 J shows 
oscillations. 

Therefore it is clearly harder to evolve adiabatically to- 
wards the expected ground state value in the presence of 
a driving force rather than in the undriven case. We at- 
tribute this finding to the many level crossings or avoided 
level crossings which occur in the quasi-energies in the 
Floquet picture. If these quasi-energy level crossings are 
too wide, the state cannot follow adiabatically and exci- 
tations are generated. For small systems it was proposed 



that a faster ramp time could, in some situations, help to 
stay in a certain state @ . Also for large systems [22| it is 
expected that short, i.e. finite, ramp time would perform 
better than completely adiabatic ramp. In our system we 
do not observe a clear signature of this phenomenon for 
the examined ramp times. In general we observe an im- 
provement of the quality of the state obtained (neglecting 
slow oscillations over the ramp time). 

Surprisingly for case (iii) the evolution of the driven 
system follows closely the evolution of the undriven sys- 
tem. This means that in this case (iii) the main devi- 
ations already stem from the internal structure of the 
energy spectrum of the undriven system. 

Additionally to the local compressibility, we investigate 
the response of the single particle correlations (jojbj+d) 
taken at the exact periods of the driving. For the ground- 
state of the superfluid (Mott-insulator) these have a char- 
acteristic asymptotic algebraic (exponential) decay with 
the distance d. In Fig. [4] (a) we show the decay of the sin- 
gle particle correlations for a fixed ramp time tf = Ylhj J. 
The target-state exhibits the exponential decay typical 
for correlations in the Mott-insulating regime. For small 
distances d the correlations decay for both the driven and 
undriven cases to a value close to the value of the tar- 
get state, whereas longer range correlations are harder to 
reach. In Fig. SJb) the dependence of the final value of 
the single particle correlations on the ramp time is shown 
for different distances d. Whereas the nearest-neighbor 
correlation (d — 1) shows a decrease towards the value 
expected for the final ground state, the correlations for 
distances d — 3 and d = 5 show an oscillating behavior 
around the target-state value. Longer ramp times would 
be necessary to converge also these correlations. 

In contrast to the compressibility, nearest-neighbor 
correlations (which are proportional to the kinetic en- 
ergy) in the driven system show very strong oscillations 
in time with a frequency uj osc = 2lu. This is exemplified 
in Fig.[5]in which the nearest neighbor correlation is plot- 
ted versus time for the case (iii), i.e. the case in which an 
effective change of the sign of J occurs. The initial value 
of the correlation is very high, since the initial state is 
close to the superfluid phase. The correlation starts to 
decrease monotonically for times corresponding to the pe- 
riod of the oscillations lo osc . However, with the increase of 
the driving amplitude strong oscillations build up. In fact 
the gas is accelerated periodically in opposite directions 
during one period of the driving potential. Only close 
to the point where the effective tunneling vanishes these 
oscillations become small in amplitude. The value at the 
exact period becomes the minimum of the oscillations for 
longer times. Due to the fast oscillations the characteris- 
tics of the state is very different if a measure is not done 
at the exact periods of the driving. While the eigenen- 
ergies of the effective Hamiltonian H e g do match closely 
the quasi-eigenenergies of the Floquet Hamiltonian fj(i), 
a lot of kinetic energy is absorbed and released during 
one period of oscillations. Hence in an experiment a very 
precise timing between the driving and the measurement 
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FIG. 4. (color online) (a) Single particle correlations between 
different sites versus their distance in the presence of the driv- 
ing amplitude of frequency Hcj/J — 25 (continuous line with 
squares) and for the undriven case (dashed line with circles). 
The target value is shown for comparison (continuous line), 
(b) Single-particle correlations (d — 1,3,5) versus the final 
time of the evolution. Parameters: U/J — 4, N = L — 48, 
j = 25 and ft such that J e a(tf) = 0.5J. Results are taken at 
exact periods of the driving. 



must be assured. 




is the time-of-flight imaging. For long time-of-flights, 
the signal of this measurement is related to the momen- 
tum distribution of the atoms which we show in Fig. [5] 
The momentum distribution is very different in the final 
target-states of the three considered situations [H,!!!)]. If 
J e ff is reduced to half its value, the initially sharp peak at 
k = broadens. At J c ff = the distribution becomes flat 
and if J c g changes sign a peak at k = ±7r/a appears. The 
final distribution for a slow ramp time is shown in Fig. [6] 
This characteristic behavior can be clearly identified in 
the momentum distribution of the driven and undriven 
evolution (taken at the exact period). However, reminis- 
cent features of the original distribution remain visible as 
a small peak at k = both in the driven and undriven sit- 
uations. The details of these features are slightly distinct 
for the driven and undriven system. In agreement with 
the deviation found for longer range correlations, the re- 
maining peak at k = is an experimentally detectable 
sign that the system has not yet reached the target-state, 
but is in an excited state. 
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FIG. 6. (color online) Momentum distribution p(k) versus 
momentum for ft such that (i) J e g(tf) = 0, (ii) J e g(t/) = 
0.5J, (iii) J e s{tf) ~ — 0.4J. The continuous line corresponds 
to the driven case, the dashed line to the undriven case and 
the stars to the value of the target-state. Parameters: U/J = 
4, hw/J = 25, N = L = 48, t f w 17H/J. 



FIG. 5. (color online) Real part of the nearest-neighbor cor- 
relation (fft((St$ i+1 ))) versus time for the driven case (dotted 
line) and undriven case (continuous line). Squares mark the 
values at each quarter of a period. Parameters: U/J = 4, 
hw/J = 25, N = L = 48 and f t such that J e s(tf) = -OA J. 

In experiments one possible measurement technique 



B. Variation of the initial state 

Until now we have focused on the case in which the 
initial state (U/J = 4) lies just beyond the critical 
value inside the Mott-insulator regime. However, since 
the superfluid and Mott-insulating state typically have 
very distinct dynamics, it is interesting to study differ- 
ent quench scenarios: from superfluid to superfluid, from 
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superfluid to Mott-insulator and from Mott-insulator to 
Mott-insulator. In Fig. [7] the final compressibility after 
a quench from J e g = J to J e g = J/2 is compared for the 
undriven, driven and target-state starting at different ini- 
tial values of U/J. The final compressibility for the un- 
driven case is always closer to the one of the target-state 
than the one of the driven case. Moreover for interme- 
diate initial values of U/J (1.7 < U/J< 3.4), that lead 
to the crossing of the phase transition, the deviations of 
both the driven and the undriven case from the target 
state are largest (see inset of Fig. We attribute 

these deviations to the change of the nature of low ly- 
ing energy levels around the phase transition. In the 
undriven case the deviations almost vanish for a quench 
inside the Mott-insulating phase and seem to be slightly 
more pronounced in the superfluid phase. The driven 
case shows larger deviations in the Mott-insulating re- 
gion than for a quench within the superfluid phase. 
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FIG. 7. (color online) Compressibility k, versus the strength 
of the interaction U. Squares represents the values of the 
driven case, circles the undriven case and the continuous line 
the value of the target-state. Inset: Deviations of the com- 
pressibility k from that of the target state kq. Parameters: 
L = N = 64, t f « 12 h/J, hid /J = 25 and f t such that 
J e g(t f ) = 0.5J. 



C. Maintaining the prepared state 

One of the purposes of using high frequency drivings 
is the preparation of complex quantum states. We have 
seen that for a slow increase of the driving amplitude 
a state close to the ground state of an effective Hamil- 
tonian can be prepared in the sense that many proper- 
ties of the driven state (measured exactly at the period 
of the driving) are close to the properties of the target- 
state. However, the question arises if this prepared many 
body state, which is continuously subject to a high fre- 
quency driving, is stable enough over a time necessary 



FIG. 8. (color online) Real part of the nearest-neighbor cor- 
relation (SR((S]S j+ i))) versus time for the driven case with 
hjj/J= 25 (dotted line, squares mark the value at a quarter 
of a period) and undriven case (continuous line). Parameters: 
U/J = 4, N = L = 48, j = 25 and f t such that (a) J eS (t) = 0, 
(b) Jeff(t) = 0.5J, (c) Jeff(t) « -0.4J for all t > tf fa 6ft/ J. 
tf is marked, in each case, with a vertical line. 



for experiments. In Fig. [8] we show the time evolution 
of the nearest-neighbor correlation for the cases in which 
the target-state is that of an effective Hamiltonian with 
J e ff = in (a), J G ff = J/2 in (b) and J e g « —0.4 J in (c). 
In all these three cases we linearly increase the amplitude 
of the driving to the desidered value in an intermediate 
ramp time tf « Qh/J. After this time the amplitude of 
the driving force is kept constant. We find that the value 
of the nearest-neighbor correlations, during this time of 
constant amplitude of the driving, oscillates around a 
quasi-steady value (when measuring exactly at the pe- 
riod of the driving). The quasi-steady state seems to be 
reasonably stable on experimentally relevant time-scales 
and no energy absorption originating from the driving is 
evident. The amplitude of the oscillations depends on the 
performed quench. In particular for the quench towards 
an effectively vanishing tunneling amplitude the oscilla- 
tions seem to be more pronounced than for the other two 
cases. These oscillations stem from the non-differentiable 
form of the protocol for the driving amplitude and are 
well known from the finite time Landau Zener transition. 
Note that the derivative of the Bessel function at its first 
zero is large, such that the abrupt change of the oscil- 
lations amplitude from linearly increasing to constant in 
time leads to the largest oscillations in case (i). The os- 
cillations might be minimized using smoother protocols 
for the change of the driving amplitude. 
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VI. CONCLUSIONS 



We have investigated for experimentally realistic sizes 
the dynamics of bosonic atoms subject to a periodic 
driving force. At small driving frequencies, the resonant 
behavior can be employed to obtain information on 
the energy spectrum of the underlying quantum phase. 
At high driving frequencies an effectively undriven 
Hamiltonian is expected to describe the properties of 
the system at the exact driving frequencies. Increasing 
slowly the driving amplitude can be used to prepare 
the corresponding ground states of the complex effective 
Hamiltonians. We have found a wide range of frequen- 
cies, depending on the properties of the underlying 
state, in which an effective Hamiltonian can be well 
simulated by the high-frequency driving. The adiabatic 
following of local quantities while increasing the driving 
amplitude, however, has been found to be less efficient 
than in the undriven case and across a phase transition. 
Longer range quantities as the single particle correla- 
tion functions are more involved to reach than local 
quantities leading to experimentally detectable signals 
in the interference patterns. Finally we have found 
that once reached the desired effective Hamiltonian, the 
created state can be maintained over a certain time 



without sizable heating. In summary, the ground-state 
of an effective Hamiltonian can be prepared to a good 
approximation using a high driving force. This is very 
promising for the engineering of more complex states by 
the application of a driving force. However, the quality 
of the prepared state is very sensitive to the preparation 
process and in particular the time employed to increase 
the amplitude of the oscillations. This will be important 
when trying to produce states with particularly low 
entropy as spin liquid states in frustrated systems. 

Note added - During the final stages of the preparation 
of this manuscript, frustrated systems have been realized 
experimentally (as show on the arXiv [III). In the light 
of these frustrated systems the question of whether 
the ground state properties of a driven system can be 
exploited becomes even more urgent. 
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